The spatial correlations in the velocities arising from a random 

distribution of point vortices. 



Pierre-Henri Chavanis and Clement Sire 

Laboratoire de Physique Quantique (UMR C5626 du CNRS), Universite Paul Sabatier 
118 route de Narbonne, 31062 Toulouse Cedex 4, France 
(chavanis@irsamc2.ups-tlse.fr & clement@irsamc2.ups-tlse.fr) 
Tel: +33-5-61558231 Fax:+33-5-61556065 
(February 1, 2008) 



Abstract 



This paper is devoted to a statistical analysis of the velocity fluctuations 
arising from a random distribution of point vortices in two-dimensional turbu- 
lence. Exact results are derived for the correlations in the velocities occurring 
at two points separated by an arbitrary distance. We find that the spatial 
correlation function decays extremely slowly with the distance. We discuss 
formal analogies with the statistics of the gravitational field in stellar systems. 
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I. INTRODUCTION 



Recently, several papers have focused on the statistics of velocity and velocity gradients 
produced by a random distribution of point vortices in two dimensions [|l]-[(J . This problem 
was first considered by Min et al. ji]] and independently by Jimenez 0. They showed, 
using different methods, that the velocity p.d.f. are Gaussian (but with a slow convergence) 
while the distribution of velocity gradients follow a Cauchy law. Their theoretical results 
were confirmed by numerical simulations of point vortices and by Direct Navier Stokes 
simulations of 2D decaying turbulence when the flow becomes dominated by a large number 
of coherent vortices 0. Additional simulations of point vortex systems were performed 
by Weiss et al. who emphasized the importance of vortex pairs in the dynamics and 
their role in the tail of the velocity p.d.f. They showed in particular that these pairs are 
responsible for anomalous diffusion and proposed to model the motion of the vortices by an 
Ornstein-Uhlenbeck stochastic process combined with Levy walks. 

A formal analogy exists between the statistics of the velocity fluctuations due to a col- 
lection of point vortices and the statistics of the gravitational field produced by a random 
distribution of stars. The only difference, of great importance, is the space dimensionality 
D = 2 instead of D = 3. In a long series of papers, Chandrasekhar and von Neumann 
|7HTTf analyzed in detail the distribution of the gravitational field arising from a random 
distribution of stars. Their study was inspired by the work of Holtsmark concerning the 
fluctuations of the electric field in a gas composed of simple ions and by the work of von 
Smoluchowski |13[ concerning the persistence of fluctuations in the Brownian theory. Chan- 
drasekhar and Von Neumann rederived the Holtsmark distribution for F, the gravitational 
field, and determined many other statistical measures for the correlations of F or the joint 
distribution of F and dF/dt. Their motivation was to derive the expression for the speed 
of fluctuations and the diffusion coefficient of stars in a purely stochastic framework. Like- 
wise in the vortex problem, the formation of binary stars alters the results of the statistical 
analysis at large field strengths. The beautiful work of Chandrasekhar and von Neumann is 
an imposing "tour de force" and provides detailed mathematical methods for analyzing the 
statistics of fluctuations in physics and astronomy. 

The close connexion between stellar and vortex systems has been investigated by one 
of us in a series of papers [TH|-[2l][] and it was natural to consider the extension of the 
Chandrasekhar-von Neumann analysis to the case of point vortices. Therefore, in Ref. 
we analyzed in some details the statistical features of the stochastic velocity field produced 
by a random distribution of point vortices in 2D turbulence. We rederived the results of 
Min et al. |l|] and Jimenez || for the distribution of velocity V and acceleration A (closely 
related to the velocity gradients) using the methods of Chandrasekhar and von Neumann. 
This formalism is very powerful and allowed us to obtain generalizations of previous works: 
we considered the distributions of the vectors V and A (not only the projection along an 
axis), we gave expressions for these distributions valid for an arbitrary number of vortices 
(for finite N these distributions can be evaluated numerically) and we extended the analysis 
to an arbitrary spectrum of circulations among the vortices and to the case of non singular 
vortex "blobs" . We also considered for the first time the joint distribution of velocity and 
acceleration ^(V, A), determined the typical duration of the velocity fluctuations T(V) 
and proposed an expression for the diffusion coefficient D of point vortices thereby justify- 
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ing the phenomenological result of Weiss et al. [[|. These theoretical results were found to 
be in good agreement with numerical simulations [21] and experiments [ 22 of 2D decaying 
turbulence. They are also of great importance to build up a rational kinetic theory of point 
vortices |T|,|2|. 

The motivation of the present paper is to characterize the spatial correlations in the 
velocities occurring at two points separated by an arbitrary distance. This problem is 
customary in turbulence but, in general, exact results are difficult to obtain. It is an interest 
of our model to allow for nice analytical solutions. Like in Ref. Q, we consider a collection of 
N point vortices with circulation 7 randomly distributed in a disk of radius R. We assume 
that the vortices have a spatial Poisson distribution, i.e. their positions are independent 
and uniformly distributed over the entire domain. We are particularly interested in the 
"thermodynamical limit" in which the number of vortices and the size of the domain go 
to infinity (N — > 00, R — > 00) in such a way that the vortex density n = remains 



finite. In this limit, the Poisson distribution is shown to be stationary [23]. Therefore, 
if the vortices are initially uniformly distributed they will remain uniformly distributed in 
average during all the subsequent evolution. This property has been checked numerically 
by Jimenez || using Direct Navier Stokes simulations of 2D decaying turbulence. In the 
statistical theory of point vortices initiated by Onsager [^3] and further developed by Joyce 
& Montgomery [[25|] and Lundgren & Pointin [2fJ, the Poisson distribution corresponds 
to a structureless equilibrium state with inverse temperature (3 = (when the angular 
momentum is not conserved) or (3 — > +00 (when the angular momentum is conserved). Of 
course, more general initial conditions are possible and lead to equilibrium states with (3 < 



in which the vortices are clustered in "macrovortices" [^,^]. The statistics of fluctuations 
remain the same in these more general situations |l],|20| but the vortices are expected to 
experience a systematic drift (JTBJ (Chavanis, 1998) in addition to their diffusive motion, 
due to the inhomogeneity of the vortex cloud. It is in this context that a kinetic theory of 
point vortices, consistent with Onsager approach, can be constructed [ IBLED ). At equilibrium 
the drift balances the scattering and maintains nontrivial density distributions. We shall 
restrict ourselves, however, in the present article to the case of a uniform distribution of 
point vortices for which the drift cancels out. 

In section |J, we recall known results concerning the distribution of velocities W(V) 
occuring at a fixed point (equations (|ICi|) (|TT|) and (|12"D). This is essentially to set the notations 
that will be used in the sequel. Since the velocity distribution created by point vortices is 
intermediate between Gaussian and Levy laws, we shall call it the "marginal Gaussian 
distribution". We also consider the case of non singular "vortex blobs" with a core of size 
a. We write down the exact characteristic function of VT(V) valid for all velocities and all 
core sizes (equation (p0[)). For "extended" vortices, we prove that the velocity distribution 
is Gaussian (with no tail) in agreement with the numerical observations of Jimenez || and 
Bracco et al. ||27|| . For a = we recover the point vortex limit for which the velocity 
distribution has a Gaussian core and an algebraic tail. For "small" non singular vortices 
(a — > but a / 0), we argue that the form of our characteristic function can explain the 
occurence of almost exponential tails observed by Jimenez and Bracco et al. |27[ . 

In section Q [II we analyze the joint distribution of velocity and velocity gradient 
W(V,5*V). In subsection [IIIB| , we rederive the Cauchy distribution for 5Y using 



the method of Chandrasekhar and von Neumann (equation (|52"D). In subsection [111 C 



we 
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generalize this method to the case of "vortex blobs" and find the explicit characteristic 
function of W(5~V) valid for all velocities and core sizes (equation (|57)). We find that the 
distribution is Cauchy for small fluctuations and Gaussian for large fluctuations. This is in 
agreement with the asymptotic behaviours found by Min et al. [p]]. It is likely that in between 
the distribution passes through an exponential tail as observed numerically in Ref. ]!]. In 
subsection LI [D], we determine an exact expression for the moment (5V)v, the average value 
of SV for a given velocity V (equations Q85p and (|5ED). When V — > +00, our results have 
a clear physical meaning in the nearest neighbor approximation. In this approximation, we 
determine the n-th conditional moments ((5V) n )v of the velocity gradients. 

In section [IV], we turn to the more difficult problem concerning the velocity correlations 
between two points separated by a finite distance. Special attention is devoted to the 
bivariate distribution W(Vq, Vi) that a velocity Vo occurs at O and that, simultaneously, 
a velocity V x occurs at a point separated from the first by a distance ri. In subsection 
IV B| , we determine an exact but complicated expression for the moment (Vi)v (equations 



m& ^d (01)). 

More explicit results are obtained in section [V] when this quantity is averaged over all 
mutual orientations of ri and Vo- An exact expression is found for (Vi)vb * n subsection 
V B| (equations ( |132| )( |i"3"4"| )). In subsection |V C[ we give the asymptotic behaviour of this 



quantity when r\ — > +00 (equations (|147|) (|HS| 

In section [V4|, we consider the projection of (Vi)v i n t ne direction of V and its average 
over all mutual orientations of ri and Vo. This defines a quantity (Vi||)y (subsection [VI A| ) 
whose expression is given by equation (|156|) . In subsection [VI B| , we give the asymptotic 
behaviours of this quantity in the limits r\ — > (equation ( |166| )) and t\ — > +00 (equations 
( |169|)(|170|) ). When (Vi||)vb i s further averaged over all possible values of V , this determines 
a function of r\ alone which characterizes the correlations in the velocities occurring simulta- 
neously at two points separated by a distance ri (subsection |VI Q ). The exact expression of 
is given by equation ( |182| ) and its asymptotic behaviours for r\ — » and n — > +00 by 
equations (|186|)(|187|) . In subsection |VI D| we calculate the correlation function (V Vi/Vq 2 ). 
Its general expression is given by equation (|191|) and its asymptotic behaviours by equations 

Finally, in section [VII] , we calculate the spatial velocity autocorrelation function and the 
energy spectrum of point vortices. We find in subsection |V11 A| that the spatial velocity 
correlation function decays extremely slowly with the distance (equation (|206|) ). The simple 
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result (V0V1) = ln(_R/ri), also derived more directly in Appendix [D], does not seem to 
have been given previously. In subsection [VII B| we use this correlation function to determine 
the energy spectrum E(k) of point vortices (equation fl213| )). This offers an alternative to 



the method considered by Novikov ||23|| . When k — > +00, we recover the classical result 
E(k) ~ k^ 1 of Novikov and for fc-^Owe find E(k) ~ k. 

Note that the spatial velocity autocorrelation function (VqVi) diverges logarithmically 
at small separations since the variance of the velocity is not defined. This is the reason 
why we consider more regular quantities such as (Vi\\) = (V0V1/V0) and (VoVx/V^ 2 ) that 
are finite at small separations. The computation of these quantities is difficult and is a 
important part of our work. 



4 



II. THE STATISTICS OF VELOCITY FLUCTUATIONS OCCURRING AT A 

FIXED POINT 



We recall in this section the results concerning the distribution of velocity fluctuations 
occuring at a fixed point. This problem was first considered by Min et al. and Jimenez 
and later on by Weiss et al. [|J, Chukbar 0, Kuvshinov & Schep || and Chavanis & Sire ||. 
Min et al. and Weiss et al. focus on the mathematical limit N — > +oo and claim that the 
velocity distribution is Gaussian by invoking a generalization of the Central Limit Theorem 
due to Ibragimov & Linnik PSl (the proof is not straightforward because the variance of the 
velocity created by a single vortex diverges logarithmically). However, they stress that the 
convergence to the Gaussian is extremely slow and that this slow convergence is responsible 
for discrepencies at large velocities. This is confirmed by direct numerical simulations of 
point vortices JTJ|§. Jimenez [Q] and Chavanis & Sire || adopt another point of view by 
treating N as a large number but not In N. This seems to be more appropriate to physical 
situations where the typical number of point vortices does not exceed 10 4 . They proved in 
that limit that the velocity p.d.f. has a Gaussian core and an algebraic tail (see equations 
(|T0[) ([TT|) and (0)). When IniV — > +oo, the tail is rejected to infinity and the distribution 



is "purely" Gaussian in agreement with the theorem of Ibragimov and Linnik. All authors 
understood that the algebraic tail is produced by the nearest neighbor and found the correct 
exponent by a phenomenological approach. A mathematical justification of this result was 
given by Jimenez || and Chavanis & Sire |J who explicitly derived the characteristic function 
(0) associated with the velocity distribution. The p 2 factor in the characteristic function 
gives rise to the Gaussian core while the term in In p gives rise to the algebraic tail @||. This 
last result is not straightforward and the rigorous proof requires to evaluate the integrals in 
the complex plane as done by Chavanis & Sire ||. The characteristic function ([?[) was also 
determined by Chukbar Q but he arbitrarily replaced the logarithmic factor appearing in 
(0) by a constant which is not justified. Therefore, he could not derive the algebraic tail. 
Kuvshinov & Schep || (2000) wrote the right distribution but justified the algebraic tail by 
the phenomenological nearest neighbor argument not by a rigorous calculation. 

Below, we recall the main results concerning the velocity p.d.f. following the presentation 
and notations of Ref. IJ. The velocity V occurring at the center O of the domain is the 
sum of the velocities 3>j (i — 1, N) produced by the N vortices: 



N 



v = 5> (i) 



i=i 

^. = SL^i (2) 

where r« denotes the position of the i th vortex and, by definition, r_i_j is the vector rotated 
by +|. Since the vortices are randomly distributed, the velocity V fluctuates. In the limit 
when the number of vortices and the size of the domain go to infinity N, R — > oo in such a 
way that the density n = remains finite, the distribution of the velocity V is given by 



W(V) = ^- 2 j A(p)e-* v d 2 p (3) 
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with 



A(p) = e~ nC W (4) 

and 

C{p) = SL[ +0 ° (l - e ^)i-ci 2 $ (5) 

Introducing polar coordinates and performing the angular integration, equation (J|) can be 
reduced to 

Because of logarithmic divergences when R — ► +oo, the "thermodynamical limit" does not 
exist. Therefore, (0) must be considered as an equivalent of Wn(V) for large N's not a true 
limit. In Ref. @||, it was found that 



For p > and iV — > oo, we have 



167T 

and for p — ► 0, we obtain 



2 2 

C(p) = ^-lniV (8) 



C(p) = -^\np (9) 



57T 



The velocity distribution (||) is just the Fourier transform of @ with (§)-©. It is found 
that (see Ref. || for detailed calculations): 

W(V)= 2 \ e -^^ v2 (V<V crit (N)) (10) 
n , -f z In TV 



where 

^rit(JV) ~ (^lniVy^ln^OniV) (12) 

As discussed in Ref. ||, the velocity p.d.f behaves in a manner which is intermediate 
between Gaussian and Levy laws. This is because the variance ($ 2 ) of the velocity created 
by a single vortex diverges logarithmically. Therefore, the Central Limit Theorem is only 
marginally applicable: the core of the velocity distribution is Gaussian while the tail decays 
algebraically as for a Levy law. Therefore, the distribution function fllPf) flTTD can be called 
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a "marginal Gaussian distribution". The velocity V cr u(N) marks the transition between the 
Gaussian core and the algebraic tail. This algebraic tail can be given a simple interpretation 
in the nearest neighbor approximation (see, e.g., Ref. ||). It is remarkable that the velocity 
created by the nearest neighbor 

V 2 (j_\ 2 j^JL nn 

n - n - ~ \2ndJ ~ 4h 2 ttR 2 1 ! 

is precisely of the same order (up to some logarithmic corrections) as the typical velocity 
created by the rest of the system 

J 7 2 \ Ar f R , . 7 2 , 2 7 2 N 



V 2 ~N(^—)~N r(r)^— d 2 r~- -IniV (14) 

WW J\ r \=d V Un 2 r 2 4nnR 2 V ; 

where d ~ n -1 / 2 denotes the interparticle distance and r(r) = \/ixR 2 the probability of 
occurrence of a vortex in r. In a sense, we can consider that the velocity is dominated 
by the contribution of the nearest neighbor and that collective effects are responsible for 
logarithmic corrections. 

The variance of the velocity can be written as 

(V 2 ) = — / -rr-^rdr (15) 
tvR 2 Jo An 2 r 2 



The integral (p~5|) diverges logarithmically at both small and large separations. This is 
consistent with the distribution flTODflTID and the formula 

r+oo 

(V 2 ) = / W(V)V 2 2nVdV (16) 
Jo 

Because of the algebraic tail, the integral (|16|) diverges logarithmically as V — ► +oo (i.e 
r — > 0). In addition, if we were to extend the Gaussian distribution fllCf ) to all velocities we 
would find that the "variance" 

<^ 2 > = ^lniV (17) 

diverges logarithmically when iV — > +oo (i.e R — ► +oo). 
The first moment of the velocity defined by 

r+oo 

(|V|) = / W(V)V27iVdV (18) 

diverges logarithmically as N — >• oo but, unlike (T^ 2 ), converges for V — ► +oo. To leading 
order in IniV, we have: 

(|V|> = (^ln^) i (19) 

If we account for a spectrum of circulations among the vortices, the previous results are 
maintained with r y 2 in place of 7 2 0. In particular, for a neutral system consisting in an 
equal number of vortices with circulation +7 and —7, the results are unchanged. This is to 
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be expected as a vortex with circulation —7 located in — r produces the same velocity as a 
vortex with circulation +7 located in r. Since the vortices are randomly distributed with 
uniform probability, the two populations are statistically equivalent. When the system is 
non neutral (but still homogeneous), the previous results derived at the center of the domain 
remain valid at any point provided that the velocity is replaced by the fluctuating velocity 
V = V — (V) where (V) = ^7vja± is the average velocity in a corresponding to a solid 



rotation |5[]. The case of inhomogeneous distribution of point vortices is considered in 

In reality, the vortices have a finite radius a which is not necessarily small (vortex 
"blobs"). If we consider that the core radius a acts as a lower cut-off, equation (^) must be 
replaced by 

c « = /f 7 ' 1 " < 20) 

This expression can be further simplified by part integrations. When a is sufficiently large, 
we can make the approximation 

and this proves that for "extended vortices" the velocity distribution is exactly Gaussian 
(with no tail) |J . This result is consistent with the numerical observations of Jimenez and 
Bracco et al. [J27J. For a = 0, we recover the marginal Gaussian distribution QTDD QTTj) with 
an algebraic tail. For "small" non singular vortices (a — > but a ^ 0), the characteristic 
function (|20|) is quadratic for p — > +00 and p — > corresponding to a Gaussian velocity 
distribution for V — > and V — > +00 (this large velocity limit is purely formal since, 
physically, the distribution must be cut at the maximum allowable velocity V max ~ 7/4-7ra 
achieved when two vortices are at distance ~ 2a from each other). However, for intermediate 
values of p, C(p) is not quadratic and depends on the value of a. This fact can probably 
explain the occurence of almost exponential tails observed by Jimenez and Bracco et al. in 
their simulations. 



III. THE VELOCITY CORRELATIONS BETWEEN TWO NEIGHBORING 

POINTS 

In this section, we analyze the spatial correlations of the velocity fluctuations. We 
first consider the correlations between two neighboring points separated by an infinitesimal 
distance Sr. 



A. The formal solution of the problem 

The difference between the velocities occurring at two points distant 5r from each other 
is given by: 

N 

5V = Y,A (22) 

i=\ 



s 



with 

n/ ( At i 9 1 y -Arlr ■ , 1 

(23) 



7 J 5r_L 2(r;5r)r; 

"27T I 



2 4 

77 rf 



where we have assumed that one of the points is at the center of the domain. This is 
similar to the expression for the acceleration A with 5r replacing v (see Ref. @, equations 
(4) (5)). The correlations in the velocities occurring between these points can be specified by 
the function W^V, 5V) which gives the simultaneous probability of the velocity V and the 
velocity increment 5V. A general expression for the bivariate probability H^(V, 5V) can be 
readily written down following Markov's method outlined in Ref. |J, section II. A. We have: 

W N (V, 5V) = J-j / A N {p,a)e-^^d 2 pd 2 a (24) 



with: 



A N (p,a) = [/\(r)e , ^r) W (25) 



and 



$=- — ^ (26) 



If we now suppose that the vortices are uniformly distributed on average, then 



r(r) = ^ (28) 



and equation fl2"5|) becomes 



^(p^)=f4^r '"^•^'^•V (29) 

\TTK z J\r\=0 / 



Since 

rR 



irR 2 J\ r \=o 

we can rewrite our expression for A^{p,a) in the form 

rlt ^ , A 



cf r = 1 (30) 



^"» = ( 1 -^/>- e, **^ 2r ) W < 31 > 



We now consider the limit when the number of vortices and the size of the domain go to 
infinity in such a way that the density remains finite 



N 

N — > oo, R — > oo, n = — — finite 

7lR Z 
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If the integral occurring in equation fl3T|) increases less rapidly than N, then 



A(p,a) 



-nC(pff) 



with 



C(p,a) 



R 



l-e l 



|r|=0 



)<fr 



(32) 



(33) 



To make progress in the evaluation of (|3"3"|), we find it more convenient to take $ as a variable 
of integration instead of r. The Jacobian of the transformation {r} — ► {$} is 



d(v) 



7 



4tt 2 $ 4 



so that 



C(p,a) 



7 



X (l_ e ^4W)) JL d 2 <f> 



(34) 



(35) 



In equation ([35]) , iff must be expressed in terms of Combining (|26"D and (|27|) we find that 

27T f _ o „ . , _1 , 

(36) 



B. The Cauchy distribution for <5V 

In this subsection, we show that the unconditional distribution of the velocity fluctuations 
<5V is governed by a 2D Cauchy law. This result was previously given by Jimenez and Min 
et al. [0] using the methods of Ibragimov & Linnik []28| and Feller p9| . On a technical point 



of view, our calculations differ from the previous authors in the sense that we consider the 
distribution of the vector <5V instead of its projection along a particular direction (the same 
remark applies to the distribution of V in section |TI|). In addition, the method presented 
here can be extended to an arbitrary spectrum of circulations among the vortices and to the 
case of vortex "blobs" (subsection |III Q ). 

According to equation (p4"|), we clearly have 



W(SV) = ~r~~~i [ A(p,a)e- i( P v ^ 5 ^d 2 pd 2 ad 2 V (37) 

l07T 4 J 

This gives the distribution of the velocity increment <5V with no condition on V. Using the 
identity 

5(x) = / e-V*d 2 p (38) 



(27T) 2 - 

the foregoing expression for W(5V) reduces to 



W(SV) = A / A(a)e~ mv d 2 a (39) 
4tt 2 J 



10 



where we have written A(p) for A(0,a). Hence, according to equations (p2|) , (|35|) and fl36|) 
we obtain 



A(a) = e 



with 



C{a) 



7 



1 



and 



atjt = -—^ a 5r ± + 2(#j_<Jr) (fcr)j 



(40) 



(41) 



(42) 



Following the usual prescription, we have let R — > oo in since the integral is convergent 
when $ — > 0. To evaluate this integral, we shall first introduce a system of coordinates with 
the x-axis in the direction of Sr. Let us denote by 6 and (3 the angles that <& and a form 
with St. Equation (|42|) now becomes 



2tt 

7 



a$ 2 |5r|sin(2#-/3) 



and we obtain 



C(a) 



With the identity 



7 



2 ,2tt 



47T 2 JO JO 



1 — COS 



— a$ 2 \Sr\sm(26) 

7 



—d<t>d6 

$•3 



2- 



cos(xsin(2#))d# = 2ttJ (j;) 



the expression fl44|) for C(<r) reduces to 



7 

With the change of variables 5r|$ 2 , we obtain 

1 r+oo 

C(a) = - 7 <r|5r| / (1 - J (x)) 

2 jo 

Integrating by parts and using the identities 

J' (x) = -Ji(x) 

and 

' + ~ JlW -<fa = 1 



x 2 



(43) 



(44) 



(45) 



(46) 



o x 



(47) 



(48) 



(49) 



we find 
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C{a) = ^\5v\a (50) 

Hence 

A(a) = e -W* (51) 

The distribution W(SV) is just the Fourier transform of (CT) . This is the 2D-Cauchy dis- 
tribution 

= 2 L |2 7 " 7sj2 ( 52 ) 



n 2 7 2 |5r| ; 

The tail of the Cauchy distribution decreases algebraically like 

msv) ~ S^i (53) 

The distribution of velocity increments W(5*V) is the same as the distribution of accel- 
erations W(A). The only difference is the occurrence of the displacement \5r\ in place of 
the average velocity |v| (see Ref. |J, equation (B30)). 

The Cauchy distribution is a particular Levy law. Accordingly, the typical velocity 
increment produced by all the vortices 

rR ^ 2 1 ^ 2 / AT \2 

SV*~Nl i r(r)f-|<5r| 2 -rf 2 r~f |5r| 2 ( — -) (54) 



/| r |=d An 2 r 4 47r V^i? 2 

is dominated by the contribution of the nearest neighbor 



0V n.n. 47r2 ^4 ^1^1 



In addition, the algebraic tail ( |55D has a clear physical interpretation in the nearest neighbor 
approximation (see, e.g., Ref. ||, section V). 

Accounting for a spectrum of circulations among the vortices and following a procedure 
similar to that adopted in the Appendix B of Ref. 0, we find that the previous results are 
maintained provided that 7 is replaced by |7|. In particular, for a neutral system consisting 
in an equal number of vortices with circulation +7 and —7, the results are unchanged. This 
is to be expected since a vortex with circulation —7 located in produces the same velocity 
increment as a vortex with circulation +7 located in r. 



C. The distribution W(5V) for vortex "blobs" 

If the vortices are not point-like but rather "blobs" , their size a acts as a lower cut-off 
and equation fliB]) must be replaced by: 

C{°) = t-l^(l- J (^v&\6r\])±dQ (56) 



2tt Jo V V 7 J J $ 3 
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After integrating by parts, we obtain 



C(a) = vra 2 
where H{x) denotes the function 



+ 2 "' l * r| "(2^' 7| * r| »' T 



H(x) = f ^-dt (58) 
Jo t 



For a — > 0, the function C(p) is quadratic: 

I07T a z 

implying that W(5V) is Gaussian for large velocity increments 

4 n ^ 4tt a 2 I r-i /■ 1 2 

W£V) ^^e"^^ 1 1 (I5VI -> +oo) (60) 

The variance of the distribution ( |5CD is 

as can be seen directly from equation ( |54l) introducing a cut-off at r ~ a instead of d. 
For cr — > +oo, the function C(cr) is linear 

C(a) ~ i 7 |5r|a (62) 

and we recover the Cauchy distribution (|52| ) for small velocity increments. Therefore the 
distribution W(8'V) makes a smooth transition from Cauchy (concave in a semi-log plot) 
for small fluctuations to Gaussian (convex in a semi-log plot) for large fluctuations. These 
asymptotic limits were found by Min et al. [|IJ, but the general expression for the charac- 
teristic function fl57|) is new. It is likely that the two distributions are connected by an 
exponential tail as predicted and observed numerically by Min et al. Of course, when a 
is reduced, the transition happens at larger fluctuations (see equation (|6l|) ) and we have 
typically a Cauchy distribution. Inversely, for "extended" vortices (large a) , the distribution 
of velocity increments tends to a Gaussian. Similar results were obtained in Ref. for the 
distribution of accelerations (note, however, that the functional form of the function C(p), 
equation (|57|), is not exactly the same). 



D. The moment (5V)v 
The average value of S~V for a given velocity V is given by 

(*V) V = / W(V, 5V)5Vd 2 (SV) (63) 
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According to equation (|24|) it can be rewritten 

W(V)(5V)v = tAz / ^(p,^)e~^ v+a5v) 5Vd 2 prf 2 ^ 2 (^V) (64) 
(27r) 4 7 

or, equivalently, 

W(V)(5V) V = ^ / ^^)^f- W }W(iV) (65) 

Integrating by parts, we obtain 

W(V)(SV) V = -i / ^(p,a) e -^™ v )rfV 2 ^ 2 (5V) (66) 

Using the identity fl38|) , we can readily carry out the integration on <5V and a to finally get: 

W(V)(SV) V = -j^- 2 1 ^{p,0)e-^d"p (67) 

To go further in the evaluation of the integral, we need to determine the behaviour of A(p,o) 
for \a\ — > 0. In Appendix ^ it is found that 

A(p,a) = e - n C(p)+iY\ Sr \^v^8)-^^m}+o(m 2 ) (68) 

where (a x , a y ) are the components of a in a system of coordinates where the x-axis coincides 
with the direction of Sr and 6 denotes the angle that p forms with Sr. We have therefore 

Mfo 0) = -C^\Sr\ sm(29)e- nC W (69) 
|^(p,0) = i^\5r\ cos(2^)e- nC ^ (70) 

L/C' y £ 

The x-component of equation fl67|) now becomes 

^(V)(TOv = —^l^rl / d0 / pc/pcos{pV A cos(x-^)}sin(2^)e- nC(rf (71) 

87T 2 JO JO 

where \ denotes the angle that V forms with Sr. Using the expansion formula ("KB) and the 



identities ( |A8[ ) (|AS|) we can carry out the angular integration to finally obtain 

W(V)(SV X ) V = T^\6r\ sin(2 X ) [ + °° J 2 {pV)e~ nC ^ pdp (72) 
Air Jo 

For V < V crit (N), the contribution of small p's in the integral ([72] ) is negligible and we can 
replace C(p) by its approximate value (§). Therefore 

W(V)(5T4) V = g|<fr| sin(2 X ) J^°° J 2 (pV)e-^ laNp2 pdp (73) 
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Using (|T0| ) and the identity 

J 2 (x)e~ ax2 xdx = 2 - 2(l + ^)e~^ (74) 



we obtain 



(5V;)v = ^(^^)|5r|sin(2 X ) (V<V cnt (N)) (75) 



2 Vn7 2 lniV 
where B(x) denotes the function 



B(x) = -(e x - 1 -x) (76) 



Similarly, we find 

n v 



{5V y ) w = -^B[£^^\5r\ C os{2 x ) (V < V crit (N)) (77) 



For V ^ V cr it(N), the integral (|72|) is dominated by small values of p and one must use 
the general expression (0) for C(p). With the change of variables z = pV, equation ( |72|) 
becomes 

W(V)(5V X ) V = £^\5r\sm(2 X ) l o + °° J2(z)e- nC ^zdz (78) 
Using the recursion formula 

Mz) = -Ji(z) - J (z) (79) 

we obtain 

W(y)(SV x )v = £^\ S r\ sin ( 2 ^{i + °° 2Ji(z)e- nC{ Vdz - £°° J (z)e~ nC ^) zd z'j (80) 

When V — > 00, the first integral is convergent while the second, equal to 2irV 2 W(V), 
decreases like V~ 2 [see equation (JTT|) ] . Therefore, to leading order in 1/V: 



W(V)(5V X ) Y ~ ^|5r| sin(2 X ) jf °° J x {z)dz (81) 
With QTTp and the identity 

/■+00 

J n (x)dx = 1 (82) 



we find 



7 

Similarly 



(5V x )v = sm(2 X ) (V > V crit (N)) (83) 
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(«yy„)v = -^^|«yr|cos(2x) {V>V crit {N)) (84) 
7 

Equations ( ff5|) (|77|) and fl33l) (Q) can be written more compactly in the form 



(5V} v = --V 2 {5r ± + 2^-^±V\ (V >V crit (N)) ( SO ) 



7 I V 2 

where _B(x) is defined by equation (|76|). Equations fl85|) flS6"D can be compared with equation 
(172) of Chandrasekhar & von Neumann |§ for the increment of the gravitational field 
between two neighboring points. 

Equation (|86|) has a clear physical interpretation in the nearest neighbor approximation. 
In Ref. g it was shown that the high velocity tail of the distribution H^(V) is produced 
solely by the nearest neighbor. Now, the velocity and the velocity difference are related to 
the position r of the nearest neighbor by [see equations (§) and 



and 



V = -^^ (87) 
2vr r 2 K ' 



2tt \ r 2 r 4 J K ' 



(5V) V = - — ^Sr ± V 2 + 2(V ± 5r)v| (89) 



Eliminating r between these two expressions, we obtain 

2ty 

7 

which coincides with the expression ([86]) of (5V)v for large V^'s. Note that 

(5V) 2 , = ^-\8r\ V 4 (90) 

T 

From fl39| ) and (p0|) we can determine all conditional moments of <5V in the nearest neighbor 
approximation (i.e., valid for large Vs). 

If we account for a spectrum of circulations among the vortices, we have in place of 
equations (pop (|3~6D : 



(5V) v = -^V 2 {5r ± + 2^^v} {V>V cnt {N)) (92) 

In particular, for a neutral system (5V)v = 0. This is to be expected (at least for a 
symmetrical distribution of circulations) since a vortex with circulation —7 located in — r 
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produces the same velocity but an opposite velocity increment as a vortex with circulation 
+7 located in r. 

Physically, we must realize that the calculations developed in this section are not valid 
for large values of |V| (for a fixed \Sr\). Indeed, large velocities are produced by a vortex 
very close to the origin. In that case, |V| = — > oo while the velocity produced in 5r is 
bound to the value |Vi| ~ 2 J s , . Therefore, the velocity difference 5V = V — Vi should 
behave like |V|. Now, according to formula (p6|), we have |<5V| ~ V 2 . This inconsistency is 
related to the fact that, when |V| — > +00, the velocity difference <5V cannot be represented 
to any degree of accuracy by the Taylor expansion (|23| ) which assumes \Sr\ <C |r|. We expect 
therefore that the results of this section will be valid only for velocities |V| < 2 -^\s T \ • m y i ew 
of these remarks, it is necessary to consider now the more general problem of the velocity 
correlations between two points separated by a finite distance. 



IV. THE CORRELATIONS IN THE VELOCITIES OCCURRING AT TWO 
POINTS SEPARATED BY A FINITE DISTANCE 



A. A general formula for W(Vq, Vi) 



The general expression for Wn(Vq, Vi), the bivariate probability that a velocity 

Vo = -fE^ (93) 
occurs at the center of the domain and that, simultaneously, a velocity 

Vi--E r4 (94) 

occurs at a point M\ separated from the first by a distance ri can be readily written down 
following Markov's method outlined in Ref. 0, section II. A. We have 

W N (V , V 1 )=*—j A N (p,a)e-^^d 2 pd 2 a (95) 

with: 

A N (p,a) = (i-—j^Ji-e l -(V+^37^) )rf 2 r j (96) 

We now consider the limit when the number of vortices and the size of the domain go to 
infinity in such a way that the density remains finite. If the integral occurring in equation 
(|96|) increases less rapidly than N, then 

A(p,a) = e - nC( fi' a) (97) 

with 

C(p,a)= (i-e^^Wyr (98) 

J\r\=0 
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Writing 

C(p,a) = C{p) + D(p,a) (99) 
where C(p) = C(p, 0) is given by equation ([?[), we have 

D(p,a) = J* (l - e~ i ^ q S E ^ L) y 2 r (100) 

Equations fl9"5] ) (P7j) ( |9"9"D and ( |100| ) formally solve the problem but it does not seem possible 



to obtain an explicit expression for D(p,a). However, like in section fD|, if we are interested 
only in the moments of Vi for a given value of V , we need only the behaviour of D(p,a) for 
\a\ — > 0. Expanding the exponential term which appears under the integral sign in a power 
series in \a\, we obtain to first order 

D(p,a) = D^(p,a) + o(\a\ 2 ) (101) 

with 

Dm ^^L^-^w dh (l02) 



According to equations ( |9~7D (|Mj) and (|101|) , A(p,o) can be expressed as: 



A(p,0) = e -^(p)-nDW(pfT) + o { pf) (1Q3) 

In Appendix [B], it is found that 

DW(p,&) = ^p(a x sm9 1 + a y cosfli) / + °° —jJ^-\ hir> 
Z7T Jo r xzirr/ 

+oo 

-(a x sin 6»x + o y cos 6>i) ^ A n (p, n) cos(n8i) 

n=l 

+oo 

-(cr^ cos 6»i - (j y svaOx) ^ B n (p,r x ) sin(n6>i) 

n=l 

7 2 7 
— — p\nR(a x sin 6*i + a y cos #i) — i—r\a y (104) 

where r> (resp. r<) is the larger (resp. smaller) of (r, ri), v4„(p, ri) and B n (p,ri) are defined 
by equations ( |B12| ) and ( |B13| ), (cr^., a y ) are the components of <r in a system of coordinates 
where the x-axis is in the direction of r x and —9i is the angle that p ± forms with 

B. The first moment (Vi)v 
The average value of Vi for a given Vo is given by: 

(V!>v = / W(V , VOV^V! (105) 
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By a procedure similar to that adopted in section [0 VP , we find that an equivalent expression 
for (Vi)vo is 



W(Vo)(V 1 ) Vo 



% r dA 



(2n) 2 J da 



{p,0)e~ i P' Vo d 2 p 



(106) 



or, using equation (|103|) 



W(Vo)(V 1 ) Vo 



in 



(2tt) 2 J da 
Introducing polar coordinates, it can be written 

in 



■ p ) e -nC(p) e -ipv 0d 2 p 



(107) 



W(V )(Vi>vo 



where 



Z(p,V )= / 
Jo 



(2tt) 2 7o 

2k Q D {1) 

~da~ 



Z(p,V )e- nC ^pdp 



(108) 



(p)e 



-ipVo cos 



(109) 



In equation (|109|) , denotes the angle that p forms with V . It is related to x an d —6i, the 
angles that Vq and p ± form with ri by: 



-X-9i 



7T 



(110) 



Therefore, an alternative expression for Z(p, Vo) is 

r2n QJ}Q) 



Z(p,V ) 
Using ( |104j ) and the identity 



o <3ff 



(p)e ipy ° sin{x+9l) ^i 



111} 



J (x) + 2 ^ J2n(x) cos(2n6) + 2% £ J 2n+1 (x) sin((2n + 1)0) (112) 



71=1 



n=0 



we find after lengthy calculations that 

r+oo ^ 

Z x (p,V ) = 17 pJi(pVo) cosx / — J: 

Jo r 



+°° dr T f 7p 

2 



27rr 



7 

lnr> — i—pJi{pVo) cosxhi-R 



+oo 



+ 7r J2\- C ^i-i(p,n) - D 2 n+i(p,r 1 )]J 2 n(pV ) sin(2nx) 

n=0 

-in^[C 2 n{p,r 1 ) - D 2 n+2(p,ri)}J2n+i(pV ) cos((2n + l)x) (113) 



n=0 



and 
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7 2 

-i—pJi{pVo) smxln-R - %-K^r 1 J Q (pV Q ) 



_7F 5Z[C2n-l(p, ri) + D 2 n+l(p, J"].)] J 2 n (pVo) COs(2nx) 
n=0 

+00 

-iirJ2[C2n(p,r 1 ) + D 2n+2 (p,r 1 )}J 2 n+i(pVo)sm((2n + l)x) (114) 



n=0 

where we have introduced the notations: 



C n (p,r 1 )=A n (p,r 1 )+B n (p,n) (115) 

D n {p,r 1 ) = A n {p,r 1 )-B n (p,r 1 ) (116) 

for n > and C n , D n = otherwise. Evaluating C n and D n as defined in equations (|115|) 
( [L16|) , we find explicitly that: 

j2 jn r+co f r <\ n j f IP 



C n (p,r 1 ) = ^p- - ^ J n+2 (117) 
Z7T n Jo r \r > / \27rr/ 



^(P.n)^r-^WS (118) 



Z7T n Jo r \r>/ \27rr 
and we recall that J- n (z) = (—l) n J n (z). 

V. THE AVERAGE VALUE OF Vi AND THE FUNCTION (Vi)y 
A. The average velocity (Vi) 
The moment (Vi) is simply obtained by averaging (Vx)v over an possible values of V : 

(Vi) = J W(V )(V 1 )v rf 2 Vo (119) 



Substituting equation ( |108| ) in equation ( |119| ) and introducing polar coordinates, this can 
be rewritten: 

(Vx) = J^yj o d X I VodVo J o pdpe- nC ^Z(p,V ,x) (120) 

where x denotes the angle that Vo forms with i^. When the integration is performed over 
X, using equations (|113|)(|114[) , it is seen that the components of (Vi) reduce to: 



(V lx ) = (121) 

and: 

20 



(Viy) = -*o / / pdpe-^inn + D 1 (p,r 1 ))J (pV ) (122) 

2 Jo Jo 

where, according to ( |118| ): 

Dl{p , ri)= - i ± p r°°j 1 (ii.y^s (123) 

27T Jo \2irr/ r> r 
Under this form, it is not possible to interchange the order of integration in (|122|). However, 



an alternative expression for (Vi y ) can be obtained along the following lines. Writing 

r+oo r+oo 

(V ly ) = dV pV J (pV )iP(p)dp (124) 
Jo Jo 



where 

n 



Hp) = -i-e- nC ^(i 1 r l + D l (p,r 1 )) (125) 
and integrating by parts with the identity 

xJo(x) = ^(xJi(x)) (126) 

we obtain: 

r+oo r+oo 

(V ly ) = - dVj P J l (pV )f(p)dp (127) 
Jo Jo 

It is now possible to interchange the order of integration. Using formula fl82|) , we find 

r+oo ^ 

(Viy) = " / ^(P)dp = V(0) - ^(+oo) = -jn (128) 
Jo I 

Hence 

(Vi> = ~n 7 ru. (129) 

This result is of course to be expected. It corresponds to the solid rotation produced by a 
uniform distribution of point vortices with density nj. 



B. The function (Vi) 



V'o 



A simple result is obtained when (V^Vq is averaged over all mutual orientations of r± 
and V . Define 

(Yi)y ° = h r (i3o) 

The integral of (Vi)v over V) is just the average velocity in ri: 

r+oo 

(V 1 )=/ H/(V )(V 1 ) yo 27T^o(i^, (131) 
Jo 
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Therefore, comparing equation ( |131| ) with equation ( |120| ) we directly obtain 

(V lx ) Vo = 



(132) 



W(V )(V ly ) Vo = -i^- I " pdpe- nC ^(i 1 r 1 + D 1 (p,r 1 ))J (pV ) (133) 



o 



Using ( |123| ), we can rewrite our expression for (Vi y )v in the form 

W(V Q )(V ly ) Vo = | o + °° p 2 dpe- nC ^J ( P V )R(p,r 1 ) (134) 

where R(p, r\) denotes the function 

7P jo \lTir J r> r 



Remembering that the range of integration in equation ( |135| ) has to be broken at r% with the 
prescription given at the end of section |IV A| , the equation defining R(p,ri) has explicitly 
the form 

R ^ ri)= -^ + Lrjj^) dr+ri pjj^L)% (13 6) 

7P r\ Jo \Zixr J J ri \lTxr J r A 

With the change of variables z = it becomes 

d/ n 27rn 2titi [2^ -ip f+°° Jxjz) 

-R(p,ri) = 1 / Ji(z)dz + - / — 5— dz (137) 

IP IP Jo 2n ri 7 J£_ z 2 

Using the identity 



d ( J n (z) 



(138) 



(139) 



dz V z n J z n 
for n = and integrating by parts, we obtain: 

IP \27rr 1 / 2nr 1 z 2 
From equation (|139|) it is readily found that 

R(p,n)^o (n->0) (140) 
i2(p,n = + — M n-+oo 141 

7P 47rrx V 7p / 
where we recall that Jifz) ~ § for 2 — > 0. 
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C. The asymptotic behaviour of (Vi)y for n — ► +00 



Substituting for the asymptotic expansion of R(p,ri) from equation (|141|) in equation 



|134j ), we obtain for r x — > +00: 

W(Vo)(V ly ) Vo = ^ [ +0 °e- nC ^Jo(pVo)pdp 

4/T JO 



The first integral is equal to 27rW / (V ) so the first term is just the average velocity {Vi y ) 
[see section |VA| . Defining 

A(Vi)y = (Vi)v - (Vi) (143) 

we have 

W(V oW V l)vo = ^ e -^ MpV ^(J£-y ip (144) 

For Vq ^ V^^iV), the integral is dominated by large values of p and we can make the 
approximation 

W(V )A(Vi>vb ^ --^^lniV f +O °e-^ laN ^J (pV )p 3 dp (145) 
64-71-5 r f Jo 



Using ( JT0| ) and the identity 



-00 



e" Qa;2 Jo(/9x)x 3 dx = -^(l-^)e-fe (146) 



2a 2 \ Aa) 



we obtain 

27T rf V n7 z m JS J 



For Vo > V cr it(N), we can follow the procedure outlined in Ref. ||, section II. C, and we 
find after some calculations 

^<Vi)v& = 7r~"~j" (Vo £ K„ t (iV)) (148) 
2vt rf 

According to these equations 

A(Vik = -f^ (149) 



A{Yi)v " = i^if {Vo ~> +oo) (150) 

This is similar to the velocity produced by a fictitious point vortex located in O with circu- 
lation +7 (when Vq — > 0) and —7 (when Vq — > +00). It is not clear whether these results 
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have a deeper significance than is apparent at first sight. We observe that A(Vi)y changes 
sign at 

v - = {^ N T (i5i) 

We can also check that 

/ ^(V )A(V42irW = (152) 
Jo 

in agreement with (|129|). 



VI. THE FIRST MOMENT OF Vi IN THE DIRECTION OF V AND ITS 

AVERAGE 

A. The moment { Vi 1 1 ) vb 

A quantity of physical interest is the average value of Vi in the direction of V . Let us 
define 

Vi\\ = Vil = V\ x cos x + V ly sin x (153) 

where lo is the unit vector in the direction of Vo. Then 

<Vi||>Vo = (Vix)v cosx+ (Vi y )v sinx (154) 

A relatively simple formula can be obtained when the foregoing expression for (Vi||}vo * s 
averaged over all possible orientations of Vo. Define 

{Vi\\)v = ^ fj {V^vM (155) 



When equation (|108|) is substituted in equation (|155|) and the integration is performed over 



X-, only the logarithmic terms survive in the series ( |113|) ( 114 ). We are thus left with 



2 

77-7 



W(V )(V Ml ) Vo = -^L I p 2 dpe- nC ^J 1 (pV )Q(p,r 1 ) (156) 



o 



where we have written 



Q(p,r 1 )= [ + ™J 2 (^-)lnr > --l\nR (157) 
Jo \2ixr J r 2 



With the usual interpretation of the notation r>, see section [IV A| > the equation defining 
Q(p, ri) has explicitly the form 

Q(p, r,) = l„r, /" J,(2L)± + / + °°^(^) tor* - ilnfl (158) 
Jo KlurJ r Jn \2iirJ r 2 

With the change of variables z = it becomes: 
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QM^J^u^ + J- 4Mi»(£)f - ilnfl (159) 

Using the identity ( |138| ) with n = 1 and integrating by parts, we obtain: 



where it is understood that e — 0. 

From equation ( |160| ), it is readily found that 

QM = l l *(^k) (ri ^ 0) (161) 

g(p,n) = -iln(^) (ri- +oo) (162) 

B. The asymptotic behaviour of (Vi||)vb for ri — > and ri — ► +oo 
First considering the limit of (K||)vb f°r r i "~ > 0, we have according to (|156|) and (|161|) 

WiVoWuU = J^e-^MpVo) In (^)p^P (163) 

With the expression (|7|) for C(p), the foregoing expression can be rewritten: 



W(V )(Vi,|>vb = ^ /_ " " C , {p)e- nG ^J l {pV )pdp (164) 



o 



Using (|126|) and integrating by parts, we find 



W(Vo)(^i||> % = ^J o +0 °e- nC ^J (pV )pdp (165) 
The integral is just 27rPU(V ). Hence 

<^i||)v„ = Vb (r x -0) (166) 

a result of course to be expected. 

Considering now the behaviour of (V\\\)v for r\ — +oo, we have according to ( |156| ) and 
(» 

W(V )(K||>vb = S ln (-) f +C °e- nC ^J 1 (pV )p 2 dp (167) 
For Vq ^ Krit(iV), we can make the approximation 

W(V )<Vi||k ^ ^ln(-) f +C °e-^^ 2 J 1 (pV )p 2 dp (168) 
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Integrating by parts, we obtain 

(Vi\\)v = i^° ln (f ) ( y o £ V crit (N)) (169) 

For Vo > V cr it(N), we can adapt the procedure outlined in Ref. section II. C and we find 
after some calculations 

(Vi\\)v = ^n(-) (V >V cnt (N)) (170) 

C. The moment (Vni) 

If we now average (Vi|[)y over Vo, we shall obtain a function of | r x | only: 

/■+oo 

(Vi||)=/ ( Vo) < 1 1 > v 2tt Vb^Vb (171) 
Jo 

which will describe the correlations in the velocities occurring simultaneously at two points 
separated by a distance |ri|. Equation ( |171| ) is clearly the same as 

(Fill) = (^p) = / ^(Vo.VOVilo^Vo^Vi (172) 
Substituting for (Vi\\)y from equation (|156|) in equation (|171|) , we have 

(Vi\\) = — tt~ / V Q dV / p 2 dpe- nC ^J 1 (pV )Q(p,r 1 ) (173) 



2tt Jo 



It is apparent that we cannot interchange the order of integration in equation (|173 ). Fol- 



lowing a method similar to that adopted in section | V A] , we rewrite equation ( |173| ) in the 
form: 



r+oo r+oo 

(Vi\\) = dV V J 1 (pV )^{p)dp (174) 
Jo Jo 



where 



Hp) = -^p 2 e- nC ^Q(p, ri ) (175) 
Integrating by parts the second integral with the identity ( f48|) we obtain 

r+oo r+oo 

(y lN > = / dV / Jo(pVo)$'(p)dp (176) 
Under this form, it is now possible to interchange the order of integration. Using formula 



32), we find that 



r+oo q>'(n) 

Vill = / -^dp (177) 
Jo p 
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Again integrating by parts we obtain 



(fill) 



o p z 



dp 



Zir Jo 



(178) 



In Appendix |C 1| , we derive an alternative form of equation (|178|) in terms of T-functions. 
Indeed, we show that 



(fill) = (|V |> 



717 



2 x i/ 2 ,+ 



.vrlnA^ 

where we have introduced the notation 



dz 



Mz) 



Z* JO 



e-H-^dt 



(7m lnA^y/2 
V 4~ 



(179) 



180) 



to measure the distance r% in terms of the average distance between vortices. Defining the 
incomplete T-function by: 



TJp+1) = r e~H p dt 
Jo 



we have 



(Vi||> = (|V„|) 



„ 7 , v<>j~°JM T ^(l Uz 



71 111 N 



(181) 



(182) 



This result can be compared with formula (117) of Chandrasekhar |10[ for the correlation 
in the gravitational forces acting at two points separated by a finite distance. 

The dependence upon r 1; through the variable s is encapsulated in the function 



/(*) 



z z 



In Appendix |C 2j , it is shown that 

f(s)~2s (s-0) 



(183) 



184) 



f(s) ~ 1.41548... + ^ Ins (s -> +oo) 

This leads to the asymptotic behaviours 

(fi||) = (iVo|>-n 7 ri + ... (n-0) 



(185) 



186) 



(fll) 



In 



|Vo|)ln(^) 



-oo 



187) 



consistent with formulae ( |166| ) and ( |1 69| ) - dTTOj ) . We observe that the velocity correlations 
decay extremely slowly with the distance r\ while the initial slope of the function (T4||) goes 
rapidly to zero at d ~ n -1 / 2 , the interparticle distance (see equation (|186|) ). 
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D. The function 

We now wish to evaluate the function 

'VoVj 



Vn 



K{ ri ) 

According to equation (|173|) , we immediately have 

Ti r ^ /* + 00 /* + 00 

K(n) = --^ dV p 2 dpe- nC{p) J 1 (pV )Q(p,r 1 ) (189) 

Z7T JO JO 

Interchanging the order of integrations and using the identity (|82|) , we obtain 

K{n) = f°° prfpe-" c ^Q(p, n) (190) 

This function is related to equation ( |178[ ) simply by introducing p in the integral. We can 
therefore repeat the steps leading to equation ( |182j ) and we obtain instead: 

K(s) = i-^J o + °°^r s ,Ai)dz (191) 

where s is defined by equation (|180|) . Explicitly, it has the form: 

w-i-ihr^ 1 -'-"* (i92 » 

It turns out that this integral can be expressed in terms of known functions. Indeed: 

Jf W = 1 -^( 1 - e "* + i ft (sr)) (193) 

where E\{z) denotes the exponential integral 

r+oo p—t 

E x {z) = / — dt (194) 

J z t 



It has the series expansion 



E 1 (z) = C - In* - £ M ( 195 ) 



n=l U[n 



where C = 0.5772... is the Euler constant. For large values of the argument z — ► +oo, we 
have: 

^)~— fl-- + 4-^ + -) (196) 



2 \ 2! Z 3 



Using these formulae, we obtain the following behaviour of K(ri) for small and large sepa- 
rations: 

K(n) = 1 - irnrl + ... (n -> 0) (197) 
K( ri ) = ^ln(£) (n-+oo) (198) 
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VII. THE SPATIAL VELOCITY CORRELATION FUNCTION AND THE 

ENERGY SPECTRUM 



A. The correlation function (Wi) 
The spatial correlation function of the velocity is defined by: 
(VoVi) = J W(V ,V 1 )V V 1 d 2 \ d 2 V 1 



(199) 



It is related to equation ( 172|) by simply introducing Vq in the integral. Therefore, according 
to equation ( |176| ), we have immediately: 



(VoVi) 



dV 



P V J (pV )^^dp 

p 



(200) 



where <3>(p) is defined by equation (|175|) . Using the identity (|126|) and integrating by parts 
the second integral, we obtain: 



(VqVO = - jf dv J Q pJl (pv Q )^-^-jdp 



dp\ p 

Inverting the order of integration and using (jS^), we find 

Jo ap\ p . 

The correlation function (|199|) therefore takes the relatively simple form: 



dp 



(VoVi) = lim 



d>'(p) 



p-o p 

Using equations ( |175|) and (|7|), the derivative of <3> is 

4vriV 



(201) 



(202) 



(203) 



$'(p) 

2 



7T 



-pe 



16tt 2 

nC(p) g(p,ri 



\n r y A p A ) 



nC(p) 



Q(p, n) 



wy^p* 

^e-^^g^p^op 2 



For p — > equation ( |160| ) reduces to: 

Q(p,r 1 ) = - 1 -\n^+0( P 2 

Therefore, according to equation ( |203| ) , we obtain: 

•< 2 'R 



(VoVaH^ln 



J: 
2tt Vri 



(204) 



(205) 



(206) 



Like the quantities ( 172 ) ( 188 ), the velocity autocorrelation function ( 206 ) decays extremely 
slowly at large distances. However, unlike ( |172j ) and (J188|), the velocity autocorrelation 



function diverges logarithmically at small separations since the variance (V 2 ) does not exist 
(see section |J). 111 Appendix [D], we give a more direct derivation of formula (|206|) . 
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B. The energy spectrum of a random distribution of point vortices 

The correlation function ( |199| ) is related to the energy spectrum of the random distribu- 
tion of point vortices by the Fourier transform 

(V(O)V(r)) = JL/ e -^^^ 2 k (207) 

Inversely, 

AltE(k) 



k 



J e ikr (V(0)V(r))d 2 r (208) 



Introducing polar coordinates and carrying out the integration over the angular variable, we 
find 

1 r+oo 

E(k) = - J o (kr)(V(0)V(r))krdr (209) 

2 jo 

According to equation (|206|) , we have for r < R: 

(V(0)V(r)> = ^ln(*) (210) 

Hence 

E(k) = — [ R J (kr) In ( -)krdr (211) 
Air Jo \r J 

Using (|126|) and integrating by parts, we get 

E(k) = ^ [ R JAkr)dr (212) 

47T JO 

Again integrating by parts with the identity (|48|) , we obtain 

2 

E(k) = ^(1 - J (kR)) (213) 
For large k the energy spectrum reduces to 

E(k) ~ ^ (k — > +oo) (214) 
which is Novikov result. At small k, using equation fl213|) , we find 



E(k) ~ ^ (fc - 0) (215) 
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VIII. CONCLUSION 



In this paper, we have analyzed in some details the statistical features of the stochastic 
velocity field produced by a random distribution of point vortices in two dimensions. In 
particular, we have obtained exact results characterizing the correlations in the velocities 
occuring at two points separated by an arbitrary distance. We have derived an explicit 
expression for the spatial velocity autocorrelation function and found that the correlations 
decay extremely slowly with the distance. The other quantities computed in this article are 
less standard quantities in turbulence but we do not see any reason why they should not 
be considered in details. They could be measured in simulations of point vortex dynamics 
and this would provide a direct confrontation between our theoretical model and a situation 
in which the temporal dynamics of the point vortices is explicitly taken into account. One 
interest of our model is to provide exact results for the velocity correlations, which is not so 
frequent in turbulence. 

Our calculations could be relevant to the context of decaying two dimensional turbulence 
when the flow becomes dominated by a large number of coherent vortices. Indeed, our model 
is based on the same assumptions as in Ref. and these assumptions have been vindi- 
cated by Direct Navier Stokes simulations and laboratory experiments. In particular, the 
assumption that the Poisson distribution is stationary is vindicated both by theoretical ar- 
guments and by the simulations of Jimenez. Of course, in decaying turbulence the density of 
vortices varies with time but it should be possible to integrate this dependance in the theory 
to make our results useful in more general situations. Also, in two dimensional turbulence 
the vortices have a finite core (vortex "blobs") and this can severely alter the predictions 
of the point vortex model. However, our formalism is general and we can introduce a lower 
cut off in the theory to take into account finite size effects. For the first time, we have 
obtained explicit expressions for the distribution of velocity and velocity increments taking 
into account the finite size of the vortices. For "extended" vortices we have proved that the 
p.d.f. of both velocity and velocity differences are Gaussian. Extended vortices occur in the 
early stage of 2D decaying turbulence (when the area covered by the vortices is still large 
H) or when the Reynolds number is low [J27[. Gaussian p.d.f. are indeed observed in these 



situations and have been discussed in Ref. [@,§7]]- At the late stages of the decline (when the 



area covered by the vortices has decreased [0) or for high Reynolds numbers [pj the point 
vortex model should be more and more accurate. It would predict a marginal Gaussian 
distribution (with an algebraic tail) for the velocity p.d.f. and a Cauchy law for the velocity 
differences. However, our analytic results show that for "small" but non singular vortices 
the p.d.f can substantially deviate from the case a = and this can possibly explain (in 
a quasi-analytical framework) the occurence of exponential tails observed and discussed in 
Ref. |[l]2"U2"7| (this point will be developed elsewhere). 

In the present study, we have exclusively considered the spatial correlations of the velocity 
occuring between two points at the same time. Another quantity of fundamental interest is 
the temporal correlation function whose integral determines the diffusion coefficient of point 
vortices through a Kubo formula p0| . It is not possible, however, to relate the spatial and 
temporal correlations functions by considering that the point vortices follow linear trajecto- 



ries with uniform velocity as is commonly done in the case of stars or electric charges [11 



Physically, the linear trajectory approximation made in plasma physics or in stellar dynam- 
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ics is not applicable here because the vortices, unlike material particles, do not have inertia 
so they do not move on their own. There is, however, a situation where such an assumption 
can be implemented. It concerns the motion of point vortices in the presence of a strong 
background shear as investigated by Chavanis |16|J20[| . In that situation we can consider, to a 
first approximation, that the point vortices follow the streamlines of the shear and evaluate 
the temporal correlation function within this assumption. The correlation function is found 
to decay as t~ 2 for t — > +00. This is a slow decay but sufficient to insure the convergence 
of the diffusion coefficient. In addition, when the vortices move in a background shear they 
experience a systematic drift |TE| normal to their mean field velocity. The drift coefficient is 
proportional to the diffusion coefficient (hence to the velocity correlation function) through 
an Einstein relation, like in Brownian theory. At equilibrium the drift balances the scatter- 
ing and maintains an inhomogeneous vortex distribution. This drift is of great importance 
to understand the organization of point vortices at negative temperatures [24 and has been 
discussed in details elsewhere |I|^D] . In the absence of background shear, we cannot eval- 
uate the temporal correlation function simply but we can estimate the diffusion coefficient 
of point vortices through the analysis of the velocity fluctuations performed by Chavanis 
& Sire 0. Therefore, these results are of great importance to build up a rational kinetic 
theory of point vortices [JTBLED . 
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APPENDIX A: DETAILED CALCULATIONS OF SUBSECTION III D 



In this Appendix, we determine the asymptotic behaviour for \a\ — > of the function 
A(p,a) defined in subsection According to equation (p2|), we need the behaviour of 
C(p,a) for \a\ — ► 0. For a = 0, C(p,o) reduces to the function C(p) introduced in section ||. 
Writing 

C(p,a)=C(p) + D( Pl a) (Al) 

we have 

^^C^ 1 - 6 ^ 2 * (A2) 

We have let — > +oo since the integral (A~2) is convergent when |<&| — > 0. For \a\ — > 0, we 



can expand the exponential term e* ^ which occurs under the integral sign in equation (|A2| 
in a power series in a. Retaining only the first term in this expansion, we have 



D(p,a) = -i» £~ e**^)^ 2 * + °(M 2 ) ( A3 ) 



Substituting for^> from equation (|3g) in equation (jA3|), we obtain 

D(p,a) = i-L / cos(^) $y ri + 2(* ± 5r)(o*) — + o(|a| 2 ) (A4) 
Ztt J\$\=o I J $ 

To evaluate this integral, we introduce a Cartesian system of coordinates where the x-axis 
is in the direction of p. We denote by (Sri, Sr 2 ) and (<7i, a 2 ) the components of 5r and a in 
this system of coordinates and we introduce 8, the angle that $ forms with p. Transforming 
to polar coordinates and setting x = p$, we obtain after some rearrangements: 

D{p,&) = i— / d6 / — cos(xcos#) 

Z7T JO JO X 

x|(o"i5r 2 + (7 2 5ri) cos(20) + (a 2 5r 2 - a x 8rx) sin(20) J + o(|a| 2 ) (A5) 
Using the expansion 

+00 

cos(xcos#) = J (x) + 2^(-l) n J 2n (:r)cos(2n#) (A6) 

n=l 

and the identities 

cos(m9) cos(n9)d9 = 7iS nm (A7) 



i-2-k 

sin(m6') sin(n^)<i6' = 7r5 nm (A8) 



'0 
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2?r 



cos(m#) sm(n9)d9 = 



(A9) 



we find that 



Since 



D(p,a) = -z 7 Mr 2 + a 2 Sn) / ^-^dx + o(\a\ 2 ) 

Jo x 



we finally obtain 



D(fi,a) 



+°° _ . , dx 

X 



•7, 



1 

2 



-t-(a 2 5r 1 + (Ti5r 2 ) + o(\a[ 



Now, according to (|3*^) ([A1|) and ( |A12| ), we have: 



A(p,a) 



-nC{p)+i^{a 2 5r 1 +a 1 5r 2 )+o(p\ 2 ) 



(A10) 



(All) 



(A12) 



(A13) 



where we recall that (5r 1 ,5r 2 ) and (01,02) are the components of 8r and a in a system of 
coordinates where the x-axis is in the direction of p. Now, to evaluate the integral ( |57| ) we 
must express 5r and a in a fixed system of coordinates independent on p. We choose this 
system such that the x-axis coincides with the direction of 5r. If 9 denotes the angle that 
p forms with Sr, the components (a x ,a y ) of a in this system are related to (01, a 2 ) by the 
transformations: 



0i = 0x cos 9 + <j v sin 9 



(A14) 



02 = 0j/ cos 9 — a x sin 9 (A15) 
Using in addition 5ri = \Sr\ cos# and 5r 2 = — \5r\ sin^, A(p,a) can be rewritten in the form 



APPENDIX B: DETAILED CALCULATIONS OF SUBSECTION IV A 

In this Appendix, we determine the expression of the function D^{p,a) defined by 
equation ( |102| ). This function can be written alternatively: 

L> (1 W) = -— l R e^T^VQn |r - ri |)rf 2 r (Bl) 
2n J\ r \=o 

Integrating by parts, we find 



D^(p,a) = ( R v(e^^ In |r - i> ± ld 2 r 
2n J\r\=o I J 



+H [ R \n\r-r 1 \a ± \/(e i ^ e ^)d 2 r (B2) 



2tt 



r =0 
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Explicitly, it has the form 



L> (1) (p,«r) = -p- I e l ^^~ In |r - Yi\a±dle r 

Z7T J|r|=B 

I* .^Inlr-r^-'^'lA (B3) 



An 2 J\ r \=o I r 2 r 4 

The first integral is a boundary term which behaves like lni? when R — > oo; we call it: 

DlLdaryM = ~ 7^ f In |r - ril^e, (B4) 



The second integral is convergent when i? — > oo; we call it: 

7 2 /-+°° c ^^i ln , ri i(gigJL 
47r 2 7| r |=o I r 2 r 4 



42^) = -^ r^tai-^ll^- 2 ^^^ (B5) 



Then 

ndV- ^ - r. 

boundary 



D«(pa) = ^L dar >,a) + Z^(p,a) (B6) 



1. The evaluation of D^ t {p,a) 

To evaluate D-^(p,a), we introduce a Cartesian system of coordinates with the a;-axis is 
the direction of p ± . We denote by (<7i, 02) the components of er in this system. Introducing 
9, the angle that r forms with p ± and transforming to polar coordinates, we obtain: 

DintM = T^P / + °°— e i2 ^ln|r-r 1 |{a 2 + 2cos^(aisin0- ( T 2 cos^} (B7) 

47T^ JO JO T 

or, equivalently, 

D$*M = t4p ^ f + °° —e l2£ ^ In |r - n|{<7i sin(20) - <r 2 cos(20)} (B8) 
47r^ Jo Jo r 

■ ~ip cos 6 

We now expand e l 2^ and ln|r — ri| which occur under the integral sign in terms of 
sinusoidal functions using the identities: 

+00 

e ixco S = Jo ^ + 2 j2 i n J n ( x ) cos(n9) (B9) 



n=l 



-00 l /r<x m 



m=l 171 ^ r > 



In|r-n| =lnr>- 53 — ( — ) cos[m(0 - 0i)] (BIO) 



where J n (x) is the Bessel function of order n, r> (resp. r<) is the larger (resp. smaller) of 
(r, ri) and $i is the angle that i"i forms with p ± . 

Using the orthogonality properties ( |A7| ) ( |A8| ) (|A9|) , we can express D^ t (p,a) in the form: 
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n (i)/ \ 7 2 [ +00 dr ( 1P \ 

+oo +oo 

-cr 2 ^ A n (p, ri) cos(n0i) - cri X! 5 ™(P> r i) sin(n6»i 



n=l 



where: 



7 2 i n /"+ 00 /r< 



„ "P 

47r n Jo 



n=l 



ry2 jn ^ +oa d r / r< \n 

Att ^ n Jo r \r 



J, 



n+2 



IP 

27rr 

7P 
1-nr 



+ 3n-2 



J, 



n-2 



IP 

f IP 
V2vrr 



(Bll) 

(B12) 
(B13) 



2. The evaluation of D 



(i) 

boundary 



Using the same system of coordinates, we can rewrite our expression for D^] undary {p,a) 
in the form: 



2t 



(1) 7 

Dboundaryifre) = ~ i ^ K R J dOe 1 "^ In |r - ri|(«7i sin0 - g 2 cos0) 



(B14) 



Substituting for (|B9|) and (|B10| ) in equation (|B14j) and carrying out the angular integration, 
we obtain 



i +oo -n /_ 



(i) / TP A 

D boundary(Pi°) = '^^{o^r) ln ^ 



7P 

27Ti? 



Jn+1 



IP 

2nR 



cos(n9i] 
sin(n0i} 



When i? — - > oo, the function -D^],^ (p, a) takes the relatively simple form: 



D boundary(P<T) = -^<7 2 ]n i2 + l~ri(-<7 2 COS X + <7i sin 01 ) 



(B15) 



(B16) 



3. Final expression of D^{p,a) 



In equations ( |B11|) and ( |B16j ), <r = (<7i, cr 2 ) is refered to a variable system of coordinates 



depending on the direction of p ± . Since W(Vq, Vj) is the Fourier transform of A(p,<r), 
related to D^(p,o) by equation ( |103| ), we need to express the vectors p and cr in a fixed 
system of coordinates. We choose this system such that the x-axis is in the direction of rj. 
The components (g x , a y ) of a in this system are related to (<Ji, g 2 ) by the transformation: 

oi = cos 0i — Gy sin 0! (B17) 

o"2 = &x sin 0i + G y cos 0i (B18) 

where it might be recalled that — 0i is the angle that p ± forms with ri. Thus, in this new 
system of coordinates, D^(p,a) has the form (|104|) . 
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APPENDIX C: DETAILED CALCULATIONS OF SUBSECTION VI C 



1. Alternative form of equation ( 178| ) 



In this Appendix, we derive an alternative form of equation ( |17<j| ) in terms of T functions. 
Substituting for Q(p,ri) from equation (|160|) in equation (|178|) , we obtain 



where 

2 



(B) = rrf r+°° o( P ) [Ok J_M 

27T J^e ^ Je Z 



dpe -nC( P ) I ^±l dz (C2) 



and we recall that e — > 0. Inverting the order of integration, we have 



or, alternatively 



/W = ^ f °° *&dz r e-^d P - ^ r dz^l f^ Z e-^d P (C4) 

7 7 

Using the identity ( |138| ), integrating by parts the first integral in equation ( P4|) and taking 
the limit e — > 0, we obtain 



r(e) = | p J^) _ i i p e _ 

2tt Uo 2 2 J Jo H 



n 'l 2 l +<X iz J J^i f^'e-^dp (05) 



2n Jo z 2 Jo 

According to equations ( |CT| ) and ( |C5| ) we therefore have 

For t*i = 0, there remains 



This is precisely the value of (|V |), see equation (|IDD, as expected from equation ( |172| ). 
Therefore, we can write 

(Vm) = (|V Q |> - ^ / + °°^^# f^'e^dp (C8) 

27T JO Z JO 

In the last integral, we can replace the function C(p) by its approximate expression @ 
without introducing any significant error. Therefore 

= <| V „|> - ^ f^ Z e-^'*d P (09) 



27T JO Z 2 







With the change of variables £ = ^ hxNp 2 , it takes the form ( |179| ). 
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2. The function f(s) 

The behaviour of f(s) for s — > can be derived from equation ( 183|) by expanding the 
incomplete T-function which occurs under the integral sign as a power series in s. To first 
order in s, we have 

rw (±) = e-H-^dt 

= jT ' (l - t + ^ + ..}jt- 1/2 dt = 2sz + o(s 3 ) (CIO) 

Thus 

f(s) = f°° ^#r s2z2 (\)dz = 2s p dz^ + ... (Gil) 
Jo z \2J Jo z 

With the identity fl49|) , we obtain 

f(s) ^2s (s - 0) (C12) 
To obtain the behaviour of f(s) for s — ► +oo, we write equation ( 183 ) in the form 

f( s ) = [ +Q °dz^- [ S2z2 e-H-^dt (C13) 
Jo z l Jo 

and take the derivative with respect to s: 

f'(s) = 2 / + °° ■h^Le-s'^dz (C14) 

JO z 

Replacing Ji(z) in the foregoing equation by its series expansion 

1 +°° (—l) n ( z\ 2n 

*<')=?SiM(i) (C15) 

and inverting the order of the integration and the summation, we obtain: 

+oo r +oo / 2 \2n 

™ = £«kTwd (0 < C16 > 

If we now introduce the variable t = s 2 z 2 , equation ( C16|) becomes 

+00 (_-\)n 1 r +oo , 



^ n!(n + l)!(2s) 2 " +1 7o 

or, alternatively 



+oo (I 1 )" 1 / 1\ 



f n!(n + 1)! (2s) 2n+1 V 2 
Integrating this series term by term, we get 
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/( s ) =i ^ + ^i ns _!f ( X )" - 1-J—r (n + (C19) 
V; 2 4^n!(n + l)!n(2s) 2 " V 2/ V ' 

where K is a constant. Numerically, we find K = 1.41548.... Therefore, to leading order in 
s, we can write 



f(s) = 1.41548... + ^- Ins (s -> +oo) (C20) 



APPENDIX D: AN ALTERNATIVE DERIVATION OF EQUATION (|20f 

The average velocity occurring in is given by the expression 



(V ) = J W(V ,V 1 )V d 2 V d 2 V 1 (Dl) 



Using equation (|95|) and integrating by parts, we obtain 



Similarly, 



r)A 8(1 
(Vo) = "^ (0 ' 0) = m ^ ( °' 0) (D2) 



r54 8(1 
(V 1 > = -i— (0,0) = m— (0,0) (D3) 



The same method applied to the correlation function 

(VoVx) = J W(V ,V 1 )V V 1 d 2 Vod 2 V 1 (D4) 



yields 



< V ° V i> = °) = "^TT^ °) + n ^( > °) (D5) 

apatj- ap oa opoa 



Hence, we can write 



(V Vi) = (V )<Vi) +71^(0,0) (D6) 

opoa 



Since (Vo) = 0, there remains 

(v » Vi> = n S (o ' o) < D7) 

Substituting explicitly for C(p,a) from equation (|98| ) in equations ( |D3| ) and ([D7|) , we 
find 

27T j |r — ri p 
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and 



n7 2 f tj_ (r - n)_L ., 



< v » v '>^y^i7^# <ir (D9) 

By part integration, the first integral can be rewritten 

(Vi) = — — z A <f\n\r — ri\e r dl (D10) 

where z is a unit vector normal to the flow and e r is a unit vector normal to the disk 
boundary at r = R. Using the expansion ( B10|) for In |r — r x | , we obtain easily 

(Vi) = ii 7 ri± (Dll) 



which corresponds to the solid rotation ( |129|) . 

Integrating by parts the second integral, we get 

(VoVi) = ^|^ln|r-n|^M/-/ 27r5(r)ln|r- ri |rf 2 r| (D12) 
Using equation (|B10| ), we obtain 

< v ° Vi > = i>(*) (D13 » 

which coincides with formula ( |206| ) . 

Of course, the relations (|D8|) and ( P9|) can be obtained independently of the formalism 
introduced in this article. Using equations (03) and (PI), we can write 



where (X) denotes the average value / T(r 1 )...T(r N )Xd 2 r 1 ...d 2 r N associated with a decor- 
related distribution of point vortices. If we split the sum in two parts, distinguishing the 
contribution of vortex pairs from the contributions of individual vortices, we obtain: 

7 2 r_u (Tj - ri)j_ \ 7 2 /JX, v (rj - r±) ± 



v,Vl) = FiS7ra + ^^") (D15) 



. ,. rt r,- — ri 



Since all vortices are identical, we find 



+ ^ N ^ N - r >I T ^ d2r I T ^Y^ d2r (D16) 

Using r(r) = and making the approximation N(N — 1) ~ N 2 for large iV's, we get 
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(VqVx) 



4^2 



r± (r - ri). 
r 2 |r — i"i I s 



-d 2 r + 



n 2 7 2 



d 2 r 



-d 2 r 



r i 



(D17) 



The second term is just the product (V )(Vi) which is equal to zero since (V ) = 0. 
Therefore, equation ( P17 ) returns our previous expression (p9|). Accordingly, the formula 
( pi3| ) for the velocity autocorrelation function can be derived very simply, independantly 
from the general formalism introduced in the article. Note, in contrast, that (Vim) and 
(VoVi/Vq 2 ) cannot be obtained by this method. Therefore, the general theory developed in 
sections |VI C| and |VI D| is necessary to compute these more complicated quantities. 
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